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Abstract 


A  new  analytical  equation  of  state  is  offered  for  metals.  The  form  is  expressible  completely 
in  terms  of  the  maximum  lattice  vibration  frequency,  macroscopically  expressible  as  the 
characteristic  temperature,  0,  and  its  derivatives.  For  one  particularly  useful  special  case,  the 
results  have  been  found  to  match  both  cold-curve  and  shock-Hugoniot  data  for  a  variety  of 
metals.  The  equation  is  of  the  Griineisen  variety,  in  which  the  Griineisen  parameter  is  specified 
as  a  function  of  volume  only.  The  equation  begins  to  lose  accuracy  in  the  pressure  range  of  one 
to  several  megabars,  perhaps  because  of  the  breakdown  of  the  Griineisen  assumption,  whereupon 
the  Griineisen  parameter  becomes  energy  dependent. 
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1.  INTRODUCTION 


An  analytical  equation  of  state,  satisfying  the  Griineisen  assumption,  was  developed 
for  metals,  which  matches  available  cold-curve  and  Hugoniot  data  for  a  variety  of  metals  out 
to  pressures  in  the  megabar  regime.  The  model  uses  the  Universal  Cold  Curve  (UCC)  model 
of  Rose  et  al.[\]2&  its  inspiration,  rather  than  being  derived  from  first  principles.  The  UCC 
model  matches  cold-compression  data  for  a  variety  of  metallic  elements,  over  wide  ranges  of 
pressure  into  the  megabar  range.  However,  unlike  the  UCC  model,  which  only  gives  pressure 
along  the  zero-degree  isotherm  (i.e.,  the  cold  curve)  in  terms  of  a  density  variable,  the 
currently  offered  model  is  expressible  in  terms  of  higher-level  thermodynamic  functions 
(namely,  the  maximum  lattice  vibrational  frequency,  (O^,  macroscopically  expressible  as  the 
characteristic  temperature,  ©(V),  and  its  derivatives).  The  current  model  therefore  allows  the 
calculation  of  state  changes  along  any  thermodynamic  path,  not  only  that  of  isothermal,  cold 
compression.  Unlike  other  follow-on  models  to  Rose  et  al,  such  as  Vinet  et  al.  [2,  3]  and 
Baonza  et  al.  [4,  5],  the  current  model’s  thermal  component  is  coupled  to  the  lattice 
compression  component  in  a  unified  way,  precisely  because  the  lattice  compression 
component  is  expressed  in  terms  of  frequency,  not  compression.  The  general  form  of  the 
ciuxent  model  requires  the  specification  of  the  characteristic  temperature  {i.e.,  lattice 
frequency)  with  volume;  however,  one  particular  form  for  0(V)  has  been  employed,  with  very 
good  resulting  agreement  to  both  cold-curve  and  shock-Hugoniot  data  for  metals.  As  a  note 
of  warning  though,  the  use  of  the  model  for  materials  undergoing  phase  change  has  not  been 
validated,  and  is  not  therefore  encouraged. 

2.  BACKGROUND 

The  UCC  model  [1]  does  a  truly  remarkable  job  of  matching  cold-compression  data 
for  a  variety  of  solid  elements,  over  wide  ranges  of  pressure  into  the  megabar  range.  The 
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UCC  model  also  predicts  the  tensile  pressure  in  materials  under  cold  expansion  out  to  rather 
large  distentions.  The  UCC  model  expresses  the  cold  energy  as 


=  £’j[l  -  (1  +  +  0.05a^)  exp(-a)] 


(1) 


where  E,,  is  the  specific  lattice  binding  energy,  and  a  is  a  dimensionless  parameter  defining 
the  relative  spacing  of  the  lattice  atoms,  given  by 


a=y\[{VIV,f^-\)  . 


(2) 


In  eqn.  (2),  the  parameter  T|  is  called  the  anharmonicity  factor  and  V  and  Vq  l^e  current 
and  ambient  specific  volumes  of  the  material,  respectively.  The  anharmonicity  factor 
represents  a  ratio  of  lengths  on  the  atomic  scale,  but  may  be  expressed  in  terms  of 
macroscopic  parameters  as 


where  Cq  is  the  bulk  sound  speed  in  the  material  at  zero  pressure  and  absolute  temperature. 
Note  that  the  use  of  the  adiabatic,  bulk  sound  speed  is  justified,  even  though  Rose  et  al  [1] 
express  their  result  in  terms  of  the  isothermal  bulk  modulus,  since  the  zero-degree  isotherm  is 
also  an  isentrope. 
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Since  the  cold  curve  is  not  only  an  isotherm  but  also  an  isentrope,  the  cold  pressure, 
p^,  is  obtainable  simply  as  the  negative  of  the  derivative  of  the  cold  energy  ( -dEJdV).  The 

UCC  model  is  essentially  a  well-reasoned  empirical  fit  to  data  derived  from  first-principle 
models  of  cohesive  energies  expressed  as  a  function  of  lattice  spacing.  By  clever  construct. 
Rose  et  al.  [1]  developed  the  nondimensional  parameter  a,  which  effectively  scales  the  cold 
curve  for  all  solids  onto  the  same  metric.  When  applied  to  available  macroscopic  cold- 
compression  data  for  a  variety  of  materials,  the  model  provides  an  extremely  accurate 
description  of  available  data,  encompassing  pressures  out  to  one  and,  in  some  cases,  several 
megabars. 

Since  the  publication  of  the  the  UCC  model  by  Rose  et  al.  [1],  there  have  been  further 
advances  following  the  lines  of  their  work.  Vinet  et  al.  [2,  3]  successfully  introduced 
temperature  effects  into  the  context  of  the  earlier  work.  It  was  shown  that  the  form  of  Rose 
was  valid,  not  only  for  zero  temperature,  but  for  any  isotherm.  The  model  of  Vinet  thus 
became  an  isothermal  equation  of  state,  and  temperature  effects  were  accounted  for  by 
modeling  the  variation  of  bulk  modulus  with  temperature.  Though  not  explicitly  casting  their 
work  in  Griineisen  form,  they  do  so  implicitly  by  citing  that  "it  has  been  observed 
experimentally  that  above  the  Debye  temperature,  [thermal  pressure]  is  a  linear  function  of 
T."  By  adding  to  this  observation  that  above  the  Debye  temperature,  thermal  energy  is 
likewise  linear  in  T,  one  may  conclude  that  V(dpldE)y  =  T  (the  Griineisen  parameter)  is  thus 
constant  for  fixed  V.  In  both  works  of  Vinet  et  al.,  the  temperature  dependence  of  the  bulk 
modulus  was  linearly  fit  to  data  for  several  substances.  With  this  linear  fitting  form,  they 
showed  that  their  equation  of  state  could  reproduce  the  Mumaghan  equation  of  state. 

Very  recent  works  by  Baonza  et  al.  [4,  5]  have  been  able  to  advance  the  field  by 
expressing  the  compressibility  of  dense  phases  using  an  inverse  power  law  as  a  function  of 
the  divergence  pressure  along  a  certain  pseudospinodal  curve  at  the  considered  temperature. 
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To  obtain  a  general  equation  of  state  at  points  off  the  pseudospinodal,  they  too,  like  their 
predecessor  Vinet,  model  the  thermal  pressure  deviations  from  this  reference  curve  using  a 
Griineisen  law,  though  in  their  case,  it  is  one  in  which  the  Griineisen  law  is  assumed  to 
follow  a  commonly  used  form — ^namely,  (F/V)  =  fitted  constant. 

A  key  distinction,  however,  between  the  work  of  Vinet  et  al.  [2,  3],  Baonza  et 
al  [4,  5]  and  what  will  be  proposed  in  the  current  work  is  that,  in  the  works  of  Vinet  and 
Baonza,  the  function  governing  the  Griineisen  reference  curve  behavior  (an  isotherm  in  the 
case  of  Vinet,  the  zero  pseudospinodal  curve  in  the  case  of  Baonza)  and  that  governing 
temperature  (Le.,  thermal  pressure)  effects  are  functionally  independent.  There  is  no  unifying 
precept  tying  the  two  functions  together.  By  contrast,  the  current  work  will  predict 
volumetric  and  thermal  effects  with  a  single  unified  function.  Volumetric  effects  are  defined 
through  the  characteristic  temperature  function,  ©,  and  thermal  effects  (via  internal  energy,  E) 
are  defined  via  the  y  function  (related  to  the  Griineisen  function  by  y=V/r),  and  are  thus 
inextricably  linked  to  0  by  definition.  In  essence,  specification  of  the  cold  curve  in  the 
context  of  the  current  model  defines  the  thermal  behavior,  or  alternately,  specification  of  the 
Gruneisen  function  for  thermal  behavior  defines  the  cold  curve.  The  key  to  this  advance  over 
prior  work  is  modeling  the  cold  curve,  not  in  terms  of  density,  as  do  Vinet  and  Baonza,  but 
doing  so  in  terms  of  lattice  frequency,  as  did  Debye  and  Einstein. 

3.  THE  GRUNEISEN  ASSUMPTION 

The  characteristic  temperature,  0,  is  an  important  thermodynamic  parameter.  The  best 
model  existing  to  date  for  it  is  the  Debye  temperature,  The  characteristic  temperature  is 
central  to  the  question  of  how  the  constant-volume  specific  heat  varies  with  temperature,  from 
the  high  temperature  limit  of  three  times  the  universal  gas  constant  per  mole  (known  as  the 
law  of  Dulong  and  Petit  [6])  to  a  value  of  zero  at  absolute  zero  temperature.  The  derivation 
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of  specific  heat  is  based  on  quantum  statistical  mechanics  and  depends  on  the  frequency 
distribution  of  vibration  modes,  ca„  within  the  crystal  lattice.  In  1907,  Einstein  proposed  that 
all  permissable  vibrational  modes  are  simultaneously  excited  with  temperature,  and  he  derived 
a  specific  heat  expression  that  was  qualitatively  correct,  but  which  approached  zero  in  too 
rapid  of  a  manner  at  absolute  zero  temperature.  Debye  remedied  this  deficiency  in  1912  by 
proposing  a  more  realistic  frequency  distribution  of  vibrational  modes  [6].  Common  to  all  of 
these  models  is  the  notion  of  a  characteristic  temperature,  by  which  the  physical  temperature, 
T,  may  be  scaled  in  order  to  achieve  the  dimensionless  temperature  parameter,  (0/7),  through 
which  the  specific  heat  is  expressed.  In  all  of  these  models,  including  that  of  Debye,  the 
characteristic  temperature  is  directly  proportional  to  the  maximum  excited  frequency  mode. 


As  Desloge  [7]  points  out  though,  the  characteristic  temperature  is  itself  a  function  of 
volume  [i.e.,  0  =  0(V)].  The  variation  of  0  with  V  is  really  a  macroscopic  manifestation  of 
the  variation  of  vibrational  frequencies  o),-  with  V.  But  since  specific  heat  theory  indicates 
that  "the  frequency  spectrum  is  entirely  determined  by  the  limiting  frequency  ...  all 
other  oscillations  change  their  frequencies  in  the  same  ratio  [as  [8]."  Thus,  the  log- 
derivative  of  any  arbitrary  vibrational  mode  with  volume  is  the  same  as  the  log-derivative  of 
the  limiting  frequency,  namely 


rflnw.  ^  _  dlno)^ 
i/lnV  dlnV 


(4) 


This  condition  quantifies  the  Griineisen  assumption.  Not  unexpectedly,  its  value  is  known  as 
the  Griineisen  parameter,  denoted  by 
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r 


Jlnco 


max 


dlnV 


(5) 


Of  course,  it  is  more  common  to  see  this  parameter  expressed  in  terms  of  macroscopic 
thermodynamic  state  variables  as 


r  =  v(dp/dE)^  .  (6) 

The  notion  that  all  vibrational  frequencies  change  with  volume  in  the  same  ratio  as  the 
limiting  frequency  is  fundamental  to  Griineisen  theory.  The  alternative  is  that  each 
vibrational  mode  would  have  its  own  associated  Griineisen  parameter,  T,-,  and  the  aggregate 
macroscopic  effect  would  necessarily  be  an  average  over  all  of  the  individual  Fj,  weighted  in 
some  manner  by  the  fractional  energy  component  associated  with  each  mode.  Such  an 
aggregate  F  would  not  only  vary  with  volume,  but  also  with  thermal  energy,  as  the  fractional 
energy  component  associated  with  each  mode  was  redistributed  with  increasing  temperature. 

In  fact,  Grodzka  [9]  summarizes  work  that  suggests  there  is  a  significant  decrease  in  the 
Griineisen  parameter  with  thermal  energy.  Similarly,  very  high  pressure  equation-of-state 
theories,  which  include  electronic  effects,  are  clear  departures  from  Griineisen  theory. 

Since,  however,  non-Griineisen  theory  is  beyond  the  scope  of  the  current  effort,  we 
limit  the  discussion  to  materials  obeying  the  constraint  of  eqn,  (4).  Since  both  the 
characteristic  temperature  and  the  Griineisen  parameter  are  related  to  the  limiting  vibration 
frequency  of  the  lattice,  (0„^  (0  directly  proportional  to  (0„^  and  F  related  through  eqn.  [5]), 

0  and  F  are,  in  fact,  related  to  each  other  as  follows: 
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0/ 


1 


(7) 


where  the  prime  denotes  differentiation  with  respect  to  specific  volume,  and  the  state  variable 
Xj/,  defined  as  V/F,  has  been  introduced  [10,  11]  for  convenience  in  manipulating  the 
thermodynamic  equations.  The  functional  form  of  the  Griineisen  parameter  will  define  the 
quantitative  relationship  expressed  in  eqn  (7).  For  example,  if  F  is  a  constant,  independent  of 
volume,  then 


0  y  =  constant  (when  F  constant). 

If,  on  the  other  hand,  xif  is  held  constant  [i.e.,  F  =  Fq/(1+p),  where  p  =  (Vo-F)/V],  then 

=  constant  (when  y  constant). 

Finally,  if  Xjr  is  linear  in  V,  such  that  V  is  held  constant  [i.e.,  F  =  Fq/(1+Pp)],  as 
Segletes  [12]  proposed  for  reasons  of  thermodynamic  stability,  then 

0'V|;'(W)  =  constant  (when  xi/  constant). 

This  last  Gruneisen  form  is  equivalent  to  that  known  as  the  Sesame  formula  [13].  Each 
unique  expression  for  F  will  produce  a  unique  algebraic  relationship  between  F  and  0,  though 
each  will  satisfy  the  differential  relationship  of  eqn.  (7). 
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The  Griineisen  equation  of  state,  which  follows  from  the  Griineisen  assumption,  may 
be  expressed  in  a  variety  of  ways,  but  amounts  to  the  following: 

.  (8) 

In  this  equation,  the  reference  functions  are  the  pressure  and  specific  internal  energy  states 
along  a  known  reference  curve,  such  as  an  isotherm,  isentrope,  or  Hugoniot.  One  of  the 
historical  problems  with  estimating  the  Griineisen  function  has  been  that  it  was  usually 
inferred  by  using  hard  experimental  data  along  one  thermodynamic  path  (e.g.,  the  Hugoniot) 
and  comparing  it  to  predictions  of  thermodynamic  states  along  another  path  [14-16]  (usually 
the  cold  curve  predicted  by  models  of  Slater  [8],  Dugdale  and  McDonald  [17],  or  others).  In 
this  sense,  previous  nonambient  Griineisen  values  were  not  even  obtained  by  indirect 
measurement,  but  rather  were  back-calculated  by  assuming  the  validity  of  a  particular  cold- 
curve  model.  Grodzka  [9]  cleverly  displayed  the  result  of  using  this  inference  approach  for 
Griineisen  estimation,  by  showing  widely  varying  estimations  of  the  Griineisen  function  for 
aluminum,  depending  upon  whose  data  and  whose  model  was  employed.  Segletes  [10-12] 
has  investigated  the  thermodynamic  stability  of  the  Griineisen  equation  and  has  been  able  to 
place  some  constraints  upon  allowable  functional  behavior  of  the  Griineisen  function.  As 
modeling  and  measurement  improve,  such  great  uncertainty  in  the  estimation  of  the  Griineisen 
parameter  has  and  will  continue  to  diminish. 

4.  THE  MODEL 

The  currently  offered  model  is  not  derived  from  first  principles  but  was  developed 
based  on  certain  functional  relationships  observed  in  the  UCC  model  of  the  cold  curve. 

Though  able  to  effectively  match  the  zero-degree  results  of  the  UCC  model,  the  current  model 
is  not,  in  fact,  an  augmented  UCC  model.  It  is,  rather,  a  complete  equation  of  state. 
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consistent  with  the  Griineisen  constraint.  By  comparison,  the  UCC  model  comprises  an 
empirical  fit  to  the  zero-degree  isotherm,  expressed  in  terms  of  a  nondimensional  lattice 
spacing.  It  embodies  mechanical  rather  than  thermodynamic  concepts. 

Since  p^^f,  and  ‘\\r  are  functions  of  volume  alone,  the  solution  to  the  Griineisen 
equation  of  state  [eqn.  (8)]  amounts  to  determining  and  G(V),  which  satisfy 

py\r-E  =  G(V)  .  (9) 

Along  the  cold  curve,  where  p^  =  -E^,  one  obtains  the  differential  equation 

y\rEi  +  E^=  -G(V)  .  (10) 

The  complementary  solution  to  eqn.  (10),  corresponding  to  G(V)  =  0,  is 

E^  =  C®  ,  (11) 

which  may  be  verified  through  the  use  of  eqn.  (7),  where  C  is  an  integration  constant. 
Incidentally,  it  is  this  complementary  solution  that  constitutes  the  cold  curve  in  both  the 
Einstein  and  the  Debye  equations  of  state  (i.e.,  no  particular  solutions  were  offered).  But  if 
Ec  involves  terms  in  0,  it  is  not  unreasonable  to  expect  the  same  of  G(V),  A  particular 
solution  to  the  equation 
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(12) 


+y  =  © 


is  y(V)  =  -0  In©. 

It  was  the  knowledge  of  these  facts  that  resulted  in  the  inspirational  leap  of  looking  at 
the  UCC  model  (eqn.  [1])  and  surmising  that  the  a  variable  was  actually  a  -iOn©  term.  Such 
a  speculation  was  also  consistent  with  eqn.  (2),  in  that  the  natural  logarithm  may,  over  a 
limited  domain,  be  well  approximated  by  a  cube-root  relation  like  eqn.  (2).  Furthermore,  it 
was  speculated  that  the  cubic  term  of  eqn.  (1)  was  a  higher-order  correction  to  an  incorrect 
functional  form  and  that  it  would  not  be  needed  if  the  logarithm  were  the  correct  form.  Such 
a  speculation  could  only  have  originated  out  of  ignorance,  of  course,  when  it  was  realized  that 
the  UCC  model  fits  data  out  to  a  =  8,  where  the  cubic  term  carries  three  times  the  weight  of 
the  linear  term.  As  a  side  note,  it  was  later  realized  that  Vinet  et  al.  [2,  3]  dropped  the  cubic 
term  of  Rose  et  al.  from  their  formulation  as  well,  since  it  complicated  the  simplicity  of  the 
form,  and  was  not  particularly  important  in  compression.  In  their  case,  however,  the  resulting 
simplification  does  adversely  affect  the  quality  of  their  fit  to  binding  energy  in  the  tensile 
region.  As  is  seen  later,  such  is  not  the  case  here,  since  the  use  of  the  logarithmic  term  does, 
in  fact,  compensate  for  the  loss  of  the  cubic  term. 

Combining  these  hypotheses  yields 

=  £jl -(1 +/)exp(-/)]  ,  (13) 


where 


/  =  -^:in(©/©o)  . 


(14) 
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Substituting  for /gives  the  cold  energy,  and  the  negative  of  its  volume-derivative  gives  the 
cold  pressure,  as 


1  -  (1  -  is:  ln(0/©„))  (©/©„/ 


(15) 


and 


Pc 


'VJ/' 


(©/©o)^  ln(©/©o)  . 


(16) 


The  general  equation  of  state  that  results  is  obtained,  using  this  cold  curve  as  the  reference  in 
eqn  (8),  to  get 


pMf-E  =  eA  (©/©o)^  -  1  +  K{K-\)  i@l%f  ln(©/©, 


(17) 


It  can  be  seen  that  eqn  (17)  becomes  greatly  idealized  if  and  when  the  parameter  K  takes  on  a 
value  of  unity.  In  order  to  satisfy  the  ambient  sound  speed  condition,  the  constant  K  must 
take  on  a  value  of 


is:  = 


(18) 


which  incidentally  is  l/(3ro)  times  the  anharmonicity  factor,  Ti,  of  Rose  et  al.  [1]  It  is  shown 
later  that  values  of  K  for  real  metals  do,  in  fact,  cluster  in  the  vicinity  of  unity.  As  such,  the 
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deviation  of  K  from  unity  tells  something  about  the  nonidealness  of  the  material.  Eqns.  (17) 
and  (18)  constitute  the  proposed  equation  of  state,  with  eqns.  (15)  and  (16)  defining  the 
associated  cold  curve.  All  that  is  required  at  this  point  is  a  specification  of  (©/©o)  as  a 
function  of  volume. 

5.  A  USEFUL  SIMPLIFICATION 

Segletes  [10-12]  has  examined  the  thermodynamic  stability  (and  instability)  associated 
with  various  functional  forms  of  the  Griineisen  parameter.  He  proposed  a  form  that  could  be 
made  not  to  violate  the  stability  metrics  that  he  examined  at  the  time.  He  expressed  the  form, 
in  terms  of  parameter  p,  as 


r  = 


To 

1+Pp 


(19) 


The  family  of  functions  that  satisfy  this  form  may  be  simply  expressed  another  way,  in  terms 
of  the  tir  variable,  as 

\|/'  =  0  .  (20) 

Assuming  that  this  form  can  adequately  characterize  metals  over  some  relevant  domain  of 
specific  volumes,  we  may  then  employ  it  to  help  us  define  our  needed  function  (©/©o).  As 
was  shown  when  solving  eqn.  (7)  for  various  Griineisen  forms,  the  effect  of  assuming  the 
validity  of  eqn.  (20)  provides  us  with  the  needed  relationship 

-  constant  ,  (21) 
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keeping  in  mind  that  \|/  and  thus  the  exponent  are  constants.  Thus,  we  have 


0/00  = 


(22) 


where  \|f  is  linear  in  V,  as  dictated  by  the  assumption  of  eqn.  (20).  Expressing  the  proposed 
cold  curve  in  light  of  eqn.  (22)  yields 


E=E 


(i- 


1  +  rW)  ln(\i//ti;„) 


1 


(23) 


and 


Pc  = 


ln(\|//\|fo)  . 

ijf'vir 


(24) 


The  equation  of  state  becomes 


p\f-E 


(\|//\i;o)(-®^>  -  l]  -  (K  -  1)  (\|f/\|;o)('®^> 


(25) 


6.  COMPARISON  TO  THE  UCC  MODEL 

Rose  et  al.  [1]  were  able  to  express  their  relationship  for  the  cold  energy  in  terms  of 
parameter  a,  which,  in  turn,  is  expressible  in  terms  of  the  anharmonicity,  t),  and  the  relative 
volume,  V/Vq.  The  correspondence  to  the  current  model  cannot  be  made  on  a  one-to-one 
basis,  as  the  current  model  has  actually  introduced  an  additional  variable.  The  term  \|f/\|fo. 
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because  of  the  assumed  linearity  of  y,  may  be  shown  to  be 


=  l-(ry)(l-wy„)  .  (26) 

Similarly,  the  constant  multiplier  of  Klx/  in  the /parameter,  eqn.  (14),  may  be  re-expressed 
in  terms  of  t|  such  that /may  be  expressed  as 


Thus,  given  values  for  the  two  parameters  T|  and  r^Xji'',  the  variables  /  [eqn.  (27)]  and  a 
[eqn.  (2)]  may  be  implicitly  related  as  a  function  of  volume.  The  cold  energy  of  the  two 
models,  eqns.  (1)  and  (23),  may  then  be  compared  as  well.  The  rQ\|/  term  determines  the 

effective  rate  of  change  of  the  Griineisen  parameter  with  volume.  For  r^^'Vi/  =  1 ,  the 

Gruneisen  parameter  is  a  constant  value  of  Fq,  while  for  =  0 ,  the  Griineisen  function 

declines  with  increasing  compression  at  its  maximum  thermodynamically  stable  rate,  as 
determined  by  Segletes  [11],  according  to  the  function  rg/(l  +  p).  In  all  cases  of  expansion, 

the  Gruneisen  parameter  asymptotes,  at  large  expansion,  to  a  value  of  l/tj/.  The  family  of 
Gruneisen  functions  determined  by  eqn.  (20)  is  plotted  versus  specific  volume  in  Fig.  1,  with 
the  parameter  varying  from  0  to  unity  in  increments  of  0.2. 

In  order  to  preserve  the  presentation  appearance  of  the  UCC  model,  data  are  plotted  as 
{EJE^  -  1)  versus  parameter  a.  In  the  first  series  of  figures  (Figs.  2a-2g),  a  family  of  curves. 
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for  the  current  model,  is  plotted  for  a  given  anharmonicity,  t),  varying  FqAj/  from  0  to  1,  at 

increments  of  0.2.  In  the  second  series  of  figures  (Figs.  3a-3e),  the  parameter  is 

specified,  and  the  family  of  curves,  for  the  current  model,  is  generated  by  varying  the 
anharmonicity  between  integral  values  of  3  and  9,  which  are  seen  to  cover  the  full  range  of 
values  noted  by  Rose  et  al  (most  elements  fall  in  the  4  to  6  range  of  anharmonicity).  In  all 
figures,  the  UCC  model  is  plotted  as  a  single  heavy  line.  Many  observations  may  be  noted 
about  the  comparison  of  the  UCC  model  with  the  currently  proposed  model.  An  enumerated 
list  is  given  as  follows: 

1)  The  general  shape  of  the  cold  curve  for  the  two  models  is  very  similar  in  form, 
thereby  indicating  that,  in  the  current  model,  the  logarithmic  function  of  /  seems  to 
compensate  for  the  now-missing  term; 

2)  From  Figs.  2,  the  curves  generated  for  the  current  model,  corresponding  to  values 
for  Fotj/  in  the  0  to  0.4  range,  match  the  UCC  model  most  poorly,  in  both  the  tensile 

(a  >  0)  and  compressive  (a  <  0)  domains. 

3)  From  Figs.  2b-2d,  for  anharmonicity  ranges  of  4  to  6  (87%  of  elements  listed  by 
Rose  et  al.  [1]  faU  in  the  ti=3.5  to  6.5  range),  a  value  of  FoV  in  the  0.7  to  0.8  range 
matches  the  UCC  tensile  model  very  well  (see  also  Figs.  3b,  3c). 

4)  From  Fig.  3c,  a  value  of  F^y^  =  0.8  matches  the  UCC  model  in  compression  for  all 
values  of  anharmonicity. 
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5)  From  Figs.  2g  and  3e,  the  higher  the  level  of  anharmonicity,  the  greater  the  value 
of  required  to  match  the  UCC  model  results  in  tension.  The  best-fit  value  for 

approaches  unity  as  t)  increases  to  a  value  of  9. 

6)  From  Fig.  2a,  for  very  low  values  of  T|  equal  to  3,  it  is  difficult  to  pick  a  single 
value  of  rpii;'  that  matches  the  UCC  results  in  tension,  the  required  value  changing 
firom  0.8  at  small  expansions  to  roughly  0.6  at  larger  expansions. 

We  may  conclude  that  the  simplification  afforded  by  the  assumption  of  eqn.  (20)  is  justifiable 
for  studying  the  behavior  of  model  over  limited  domains  of  compression  and  expansion,  based 
on  the  quality  of  the  match  to  the  UCC  model.  The  general  model,  however,  retains  the 
flexibility  to  permit  other  functional  forms  on  the  Griineisen  parameter  and  will  be  explored 
more  fully  in  the  future. 

7.  COMPARISON  TO  DATA 

The  current  model,  incorporating  the  simplifying  assumption  of  eqn.  (20),  has  been 
used  to  smdy  compression  data  for  various  materials.  Since  the  current  model  is  a  full 
equation  of  state,  its  predictions  may  be  compared  not  only  for  cold-curve  compression,  but 
also  for  other  thermodynamic  processes  (e.g.,  shock  loading,  for  which  ample  data 
exists  [18,  19]). 

Though  it  can  be  done,  the  comparison  intentionally  avoids  taking  the  approach  of 
tuning  all  available  parameters  so  as  to  best  correlate  the  model  to  both  cold  and  shock  data 
for  several  reasons.  First,  it  avoids  the  criticism  of  so-called  'Tcnob  turning"  to  match  data. 

Of  course,  the  current  model  does  have  a  single  parameter,  and  that  is  the  value  of 
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Secondly,  it  is  known  that  the  Griineisen  assumption  becomes  invalid  at  large  enough 
pressures,  and  so  an  examination  of  when  the  current  model  diverges  from  the  shock  data 
could  be  used  to  indicate  something  about  the  range  of  validity  of  the  Griineisen  assumption. 

When  possible,  model  parameters  were  derived  directly  from  the  molecular  weights  as 
well  as  the  t|  and  hE  values  given  by  Rose  et  al.  [1].  These  derived  parameters  include  Ef, 
and  Cq.  With  the  addition  of  the  handbook  value  for  the  thermal  expansion  coefficient,  a, 
values  for  Fq  and  K  are  directly  obtained.  Densities  were  set  to  the  values  from  which  most 
or  all  of  the  experimental  data  were  obtained. 

The  only  parameter  that  was  subject  to  fitting  was  the  value  of  FoV,  which 

determines  the  manner  in  which  the  Griineisen  function  changes  with  volume  (recall  Fig.  1). 
This  parameter  was  fitted  not  to  give  the  best  overall  match  to  both  the  shock  and  cold-curve 
data.  Rather,  this  parameter  was  fitted  to  give  the  best  fit  to  the  cold-curve  alone.  In  this 
way,  the  model’s  Hugoniot  is  merely  a  prediction  of  and  not  an  explicit  fit  to  the  Hugoniot 
data.  This  distinction  is  important,  as  better  fits  to  the  aggregate  cold/Hugoniot  data  are 
obtainable  by  varying  values  of  Q  and  Fq  within  previously  published  ranges. 

The  model  data  are  collected  in  Table  1.  The  parameter  W  is  simply  the  molecular 
weight,  and  is  used  for  converting  molal  quantities  to  a  specific  (i.e.,  per-mass)  basis.  The 

second  to  last  column,  Foil;',  is  the  only  free  parameter.  Its  value  reflects  the  best  fit  to  the 

cold-curve  data.  In  the  last  column,  \)r',  which  is  also  dimensionless,  is  derived  from  other 
column  data  and  is  provided  for  informational  purposes  only.  Graphical  results  are  shown  in 
Figs.  4-12  for  a  variety  of  metals.  The  lower  curve  is  the  fitted  cold  pressure  curve.  The 
upper  curve  (when  present)  is  the  predicted  Hugoniot  curve.  Cold  data  (including  theoretical 
predictions  thereof)  depicted  in  the  figures  came  from  Mao  et  al  [22],  or  were  digitized  from 
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the  graphs  of  Rose  et  ah  [1]  and  Molodets  [23].  All  shock  data,  against  which  the  model 
predictions  were  compared,  were  taken  from  the  Livermore  compendium  [18],  the  Kohn 
compilation  [19],  or  Mitchell  et  al.  [24].  All  cold  data  are  depicted  with  solid-filled  circular 
symbols.  All  other  (unfilled)  symbols  represent  Hugoniot  data. 


Table  1.  Model  Parameters 


■ 

Pc 

(kg/m’) 

w 

(g/mol) 

'n 

E 

(MJ/kg) 

C  ^ 

(m/s) 

a*  (xlO*) 
(°C‘) 

F 
^  0 

K» 

ToV 

y 

Ag 

10490 

107.868 

5.94 

2.647 

3221 

19.0 

2.55 

0.78 

0.84 

0.33 

Al 

2700 

26.982 

4.71 

11.939 

5424 

25.0 

2.39 

0.66 

0.81 

0.34 

Cu 

8930 

63.546 

5.20 

5.312 

3995 

16.6 

2.02 

0.86 

0.78 

0.39 

Li 

537 

6.941 

3.10 

22.928 

4948 

- 

- 

- 

0.81 

- 

Mo 

10200 

95.94 

5.85 

6.846 

5102 

5.0 

1.50 

1.30 

0.43 

0.29 

Pd 

12020 

106.4 

6.41 

3.568 

4036 

11.7 

2.44 

0.88 

0.61 

0.25 

Pt 

21430 

195.09 

6.47 

2.893 

3515 

9.0 

■  2.61 

0.79 

0.78 

0.30 

Rb 

1530 

85.47 

4.18 

0.968 

1371 

- 

• 

- 

0.88 

- 

B 

7896 

- - 

56.05 

5.24 

7.342 

4733 

12.0 

1.81 

0.97 

0.76 

0.42 

^Ei,  calculated  from  AE  values  of  Rose  et  al.  [1]  using  Ej(MJ/kg)  =  96.45  A£(eV)/W. 

®Co  calculated  using  eqn.  (3)  and  values  of  n  given  by  Rose  et  al.  [1],  rather  than  handbook  values.  Platinum  is  the 

exception,  where  the  value  of  3668  m/s  yielded  by  eqn.  (3)  produces  a  fit  to  the  cold  curve  that  is  not  as  good  as 
when  using  the  handbook  value  of  3515  m/s. 

^Values  for  thermal  expansion  coefficient  taken  from  Weast  [20],  except  for  palladium,  which  was  taken  from  Layne  [21], 
and  lithium  and  rubidium,  which  could  not  be  located  in  the  literature. 

^^For  consistency,  initial  Griineisen  values  computed  using  Fq  =  aC^WIR,  where  R  is  the  universal  constant,  equal  to 

8314  g  mVs^  mol  K.  The  values  for  Fq  are  generally  compatible  with  published  values  found  elsewhere.  Values  of 
a  for  lithium  and  rubidium  could  not  be  located  (nor  could  estimates  of  Fq). 
calculated  using  eqn.  (18). 

’^^operties  for  304  stainless  steel  calculated,  per  Rose  et  aL  [1],  using  70%  Fe,  19%  Cr,  ll%Ni. 
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The  fits  to  the  cold-compression  data  are  extremely  good,  with  the  model  curve 
literally  overlaying  the  data  in  nearly  every  case.  For  most  cases,  the  Hugoniot  predictions 
are  very  good  to  a  point,  though  model  results  in  some  cases,  notably  silver  and  platinum, 
diverge  from  data  at  higher  pressures.  Lithium  and  rubidium  are  special  cases,  in  that  it  does 
not  appear  possible  that  the  cold  [1]  and  Hugoniot  [18]  data  were  even  generated  from  the 
same  density  base  material.  In  particular,  Hugoniot  data,  which  should  universally  lie  at 
higher  pressures  than  the  cold  curve,  are  not  seen  to  exceed  the  cold-curve  data  until  a 
relative  compression  of  approximately  2.  For  this  reason,  and  the  fact  that  thermal  expansion 
data  and  Griineisen  estimates  could  not  be  obtained  for  the  materials,  no  attempt  was  even 
made  to  model  the  lithium  and  rubidium  Hugoniots.  In  the  case  of  aluminum,  the  poor  fit  at 
even  moderate  pressures  is  perhaps  attributable  to  a  phase  change  since,  as  shown  in  Figure  5 
with  the  dotted  curve,  when  using  an  initial  Gruneisen  value  of  1.2  (half  of  the  value  shown 
in  Table  1),  an  excellent  fit  to  the  aluminum  Hugoniot  prevails  to  9  Mbar.  Such  a  notion  is 
buttressed  by  the  work  of  McKenna  and  Pastine  [16],  whose  empirical  fit  to  the  aluminum 
Gruneisen  function  decreases  with  compression  at  such  a  prodigious  rate  (in  the  V/Vq  range 
from  1.0  to  0.85)  that  it  violates  the  stability  metrics  of  Segletes  [10,  11],  established  for 
thermodynamic  transitions  of  single-phase  Gruneisen  materials.  Thus,  if  McKenna’s  fit  is 
accurate,  then  either  the  aluminum  Gruneisen  function  varies  with  temperature  (violating  the 
Gruneisen  assumption)  or  else  the  aluminum  Hugoniot  experiences  a  phase  change  in  VIVq 
range  of  0.85  to  1.0. 

Table  2  gives  a  very  approximate  estimation  of  the  goodness  of  the  match  of  the 
current  model  to  the  Hugoniot  data,  in  terms  of  pressure,  and  includes,  in  the  last  column,  the 
Debye  temperature  for  each  material.  In  all  cases  except  aluminum,  good  matches  to  the 
Hugoniot  are  obtained  out  to  pressures  of  approximately  1  Mbar  or  greater.  Furthermore,  the 
two  materials,  silver  and  platinum,  that  were  observed  to  experience  only  fair  correlations  at 
higher  pressures,  happen  to  be  the  materials  with  notably  low  Debye  temperatures.  It  is 
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Table  2.  Limiting  Pressure  Ranges  (Mbar)  to  which  Model  Prediction  Qualitatively 
Correlates  with  Hugoniot  Data 


Quality  of  Model  Correlation  with  Hugoniot  Data 

©o' 

Excellent 

Good 

Fair 

(K) 

Ag 

0.5 

0.9 

1 .6+ 

220 

A1 

— 

0.5 

0.7 

380 

Cu 

1.2 

4.2+ 

— 

310 

Mo 

0.6 

1.7+ 

— 

375 

Pd 

0.5 

2.2+ 

— 

— 

Pt 

2.0 

2.0 

2.7+ 

225 

304  St. 

Steel 

1.9+ 

— 

— 

433* 

^ata  extracted  from  Andrews  [6]. 
*Value  reflects  that  of  pure  iron. 


surmised  that  perhaps  the  breakdown  of  the  Griineisen  assumption,  eqn.  (4),  at  high  pressures 
(actually  temperature)  is  responsible  for  the  divergence  of  the  model  from  the  data.  In  the 
case  of  aluminum,  a  possible  cause  of  the  poor  Hugoniot  prediction  has  already  been 
discussed.  In  line  with  specific  heat  theory  and  the  notion  of  a  characteristic  temperature,  it 
would  seem  to  make  proper  thermodynamic  sense  that  the  breakdown  of  the  Griineisen 
assumption  with  temperature  might  also  scale  with  the  characteristic  temperature  [i.e.,  the 
breakdown  occurs  as  a  function  of  (0/7)].  With  the  exception  of  aluminum,  the  results  of 
Table  2  would  seem  to  support  this  reasoning. 
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8.  CONSTRAINTS  ON  THE  FREE  PARAMETER,  Toll/ 


The  free  parameter  in  the  current  model,  Tg  V,  obviously  has  physical  meaning  and 

relates  to  the  rate  of  Griineisen  parameter  change  with  volume.  Ideally,  it  could  be  quantified 
from  other  physical  constants.  Since  the  characteristic  temperature  in  the  current  model  is 
always  normalized  by  its  ambient  value  though,  one  is  limited  in  trying  to  extract  additional 
information  directly  from  this  form.  The  assumed  Griineisen  form,  eqn.  (20),  dictates  that 

=  1  be  an  upper  limit  on  the  model  parameter.  For  a  lower  limit,  stability  analyses 

conducted  by  Segletes  [11]  showed,  on  the  basis  of  high  pressure  isobaric  transitions,  that  the 
maximum  rate  of  descent  of  the  Griineisen  function  with  compression  is  limited  such  that 

\|;'  >  0.  At  cold  pressures,  however,  a  different  isobaric  criterion  arises,  in  which  it  was 
shown  that  a  stable  compressive  cold  curve  must  obey  the  relation 

-p'/p^>l/\|r  .  (28) 

Applying  this  criterion  to  the  current  model  gives  a  lower  bound  on  the  model  parameter  as 


\(/  > 


1  -K 


ln(t|r/\iro) 


(29) 


It  should  be  immediately  noted  that  the  value  of  the  parameter  K,  which  equals  ti/SFo,  drives 
the  lower  bound  on  \(r'.  In  particular,  for  materials  in  which  K  is  greater  than  unity,  the  high- 
pressure  'V}/  >  0  criterion  is  the  most  restrictive  constraint.  For  K  less  than  unity,  the  low- 
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pressure  constraint,  given  by  eqn.  (29),  governs  the  problem.  For  K  equal  to  unity,  both 
criteria  coincide  and  the  equation  of  state,  eqn.  (17),  simplifies  immensely,  with  the  right- 
hand  term  identically  zero  and  all  exponents  equal  to  unity. 

Considering  eqn.  (29),  for  the  worst-case  compression  scenerio,  one  may  take  the 
argument  of  the  logarithm  as  (1  -  r^ti/).  The  results  of  this  criterion  give  lower  allowable 

limits  on  the  model  parameter,  below  the  fitted  value  for  each  of  the  materials  studied, 

with  minimum  levels  generally  below  the  0.4  level.  This  result  reinforces  our  comparison 
with  the  UCC  model,  in  which  low  values  (0  to  0.4)  of  the  model  parameter  did  not  fit  the 
UCC  model  in  tension  well.  The  span  from  0.4  to  1,  though,  still  leaves  a  large,  permissable 

range  on  the  value  of  the  parameter  Further  analysis  may  serve  to  better  define  this 

parameter,  either  in  terms  of  more  restrictive  limits  or  in  terms  of  fundamental  quantity 
groupings. 

9.  CONCLUSIONS 

An  equation  of  state  for  metals  has  been  proposed.  The  form  of  the  equation  provides 
a  very  close  match  to  the  cold-temperature  expansion  models,  cold-compression  data,  and 
shock-compression  data.  The  model  was  inspired  by  the  Universal  Cold  Curve  (UCC)  model 
of  Rose  et  al.  [1],  when  it  was  realized  that  algebraic  terms  in  the  UCC  model  bore  a 
functional  resemblence  to  thermodynamic  terms  that  could  be  the  solution  to  the  Griineisen 
cold  curve.  In  the  model,  the  cold  energy  is  completely  expressible  with  terms  involving  the 
material’s  characteristic  temperature. 

A  simplified  form  of  the  current  model  has  one  free  parameter  per  given  material, 
which  physically  relates  to  the  rate  at  which  the  Griineisen  parameter  changes  with  volume. 
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The  fitting  parameter  for  the  model  is  not  totally  free,  but  is  constrained  within  certain 
boimds,  based  on  considerations  of  thermodynamic  stability.  With  this  one  parameter, 
excellent  fits  to  cold  curves  are  produced  to  all  available  low-temperature  compression  data, 
including  platinum  cold  compressions  to  6  Mbar.  Once  the  fitting  parameter  is  set  for  the 
cold  curve,  other  thermodynamic  paths  may  be  predicted  without  the  need  of  additional  fitting 
parameters.  The  resulting  predictions  of  Hugoniot  behavior  match  shock-compression  data 
well,  typically  to  the  megabar  range,  before  diverging.  The  divergence  of  the  current  model 
with  very  high  pressure  shock  data  might  be  attributable  to  the  breakdown  of  the  Griineisen 
assumption  at  high  pressures. 

Additional  avenues  remain  to  be  explored  with  the  current  model,  which  include 
enhancing  the  formulation  to  properly  account  for  the  very  high  pressure  regions,  examining 
more  fully  the  properties  of  the  equation  of  state  in  tension,  and  examining  the  general  form 
of  the  equation  of  state,  eqn.  (15).  Ultimately,  one  seeks  to  understand  the  underlying  physics 
of  lattice  compression  and  vibration  that  would  lead  to  an  equation  of  state  of  the  form 
proposed.  As  it  is,  the  current  model  hopefully  serves  to  put  the  often-used  Griineisen 
equation  of  state  on  a  firmer  foundation. 
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a 

(a) 

Figure  2.  Comparison  of  binding  energy  of  current  model  to  that  of  Rose  et  al.  [1],  for  given 
anharmonicities.  Rose  curve  depicted  with  dark  line.  Current  model  shown  in  dashed  lines 
for  incremental  values  of  ,  ranging  from  0  to  1  in  increments  of  0.2.  (a)  T\  =  3. 

(b)  r\  =  4.  (c)r]=  5.  (d)  ^  =  6.  (e)r]  =  7.  (0^]  =  8.  (g)r]  =  9. 
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(b) 
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(c) 

Figure  2  (continued).  Comparison  of  binding  energy  of  current  model  to  that  of  Rose 
et  al  [1  ],  for  given  anharmonicities.  Rose  curve  depicted  with  dark  line.  Current  model 

shown  in  dashed  lines  for  incremental  values  of  FgH/,  ranging  from  0  to  1  in  increments 
of  0.2.  (a)  Ti  =  5.  (b)x\=  4.  (c)  Tj  =  5.  (d)r]  =  6.  (e)r\  =  7.  (/)  t)  =  5.  (g)  T|  =  9. 
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(e) 

Figure  2  (continued).  Comparison  of  binding  energy  of  current  model  to  that  of  Rose 
et  al.  [1],  for  given  anharmonicities.  Rose  curve  depicted  with  dark  line.  Current  model 

shown  in  dashed  lines  for  incremental  values  of  rQ\|/',  ranging  from  0  to  1  in  increments 
of  0.2.  (a)  Ti  =  5.  (b)r\  =  4.  (c)  r\  =  5.  (d)  t]  =  6.  (e)r]  =  7.  (/)  t|  =  8.  (g)  T)  =  9. 
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(8) 

F^ure  2  (continued).  Comparison  of  binding  energy  of  current  model  to  that  of  Rose 
et  al.  [1  ],  for  given  anharmonicities.  Rose  curve  depicted  with  dark  line.  Current  model 

shown  in  dashed  lines  for  incremental  values  of  Tp-vif',  ranging  from  0  to  1  in  increments 
of  0.2.  (a)  11  =  5.  (b)T]  =  4.  (c)  r]  =  5.  (d)r\  =  6.  (e)y\  =  7.  (f)r\  =  8.  (g)  r]  =  9. 
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(b) 

Figure  3.  Comparison  of  binding  energy  of  the  current  model  to  that  of  Rose  et  al.  [1],  for 
given  values  of  Fq  \|/  parameter.  Rose  curve  depicted  with  dark  line.  Current  model  shown 
in  dashed  lines  for  integral  values  ofr\,  ranging  from  3  to  9. 

(a)  =  0.6.  (b)  =  0.7.  (c)  =  O.S.  (d)  =  0.9.  (e)  F^V  =  1.0. 
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a 

(d) 

F^ure  3  (continued).  Comparison  of  binding  energy  of  the  current  model  to  that  of 
Rose  et  al.  [1],  for  given  values  of  parameter.  Rose  curve  depicted  with  dark  line. 
Current  model  shown  in  dashed  lines  for  integral  values  ofr\,  ranging  from  3  to  9. 

(a)  =  0.6.  (b)  =  0.7.  (c)  =  0.8.  (d)  V  =  0.9.  (e)  =  1.0. 
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(e) 

Figure  3  (continued).  Comparison  of  binding  energy  of  the  current  model  to  that  of 
Rose  et  al.  [1],  for  given  values  of  parameter.  Rose  curve  depicted  with  dark  line. 
Current  model  shown  in  dashed  lines  for  integral  values  o/r),  ranging  from  3  to  9. 

(a)  =  0.6.  (b)  ToV  =  0.7.  (c)  V  =  0.8.  (d)  F^V  =  0.9.  (e)  F^V  =  1.0. 


va/q 

Figure  4.  Cold-curve  fit  with  Hugoniot  prediction  for  silver,  to  2  Mbar.  Cold  data  [22] 
(filled  symbols)  and  Hugoniot  data  [18, 19]  (unfilled  symbols)  included. 
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Figure  5.  Cold-curve  fit  with  Hugoniot  prediction  for  aluminum,  to  9  Mbar.  Cold  data  [23] 
(filled  symbols)  and  Hugoniot  data  [19,  24]  (unfilled  symbols)  included.  Dashed  curve  is 
Hugoniot  prediction,  assuming  Fq  =  1.2. 


V/Vo 

Figure  6.  Cold-curve  fit  with  Hugoniot  prediction  for  copper,  to  5  Mbar.  Cold  data  [22] 
(filled  symbols)  and  Hugoniot  data  [19]  (unfilled  symbols)  included. 
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Figure  7.  Cold-curve  fit  for  lithium,  to  0.8  Mbar.  Cold  data  [1]  (filled  symbols)  and 
Hugoniot  data  [18]  (unfilled  symbols)  included. 


MN, 

Figure  8.  Cold-curve  fit  with  Hugoniot  prediction  for  molybdenum,  to  2  Mbar.  Cold 
data  [22]  (filled  symbols)  and  Hugoniot  data  [19]  (unfilled  symbols)  included. 
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Figure  9.  Cold-curve  fit  with  Hugoniot  prediction  for  palladium,  to  2.5  Mbar.  Cold 
data  [22]  (filled  symbols)  and  Hugoniot  data  [18]  (unfilled  symbols)  included. 


Figure  10.  Cold-curve  fit  with  Hugoniot  prediction  for  platinum,  to  8  Mbar.  Cold  data  [1] 
(filled  symbols)  and  Hugoniot  data  [19]  (unfilled  symbols)  included. 
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Figure  11.  Cold-curve  fit  for  rubidium,  to  0.4  Mbar.  Cold  data  [1]  (filled  symbols)  and 
Hugoniot  data  [18]  (unfilled  symbols)  included. 


.5  .6  .7  .8  .9  1.0 


V/Vo 

Figure  12.  Cold-curve  fit  with  Hugoniot  prediction  for  304  stainless  steel  (70%  Fe,  19%  Cr, 
11%  Ni),  to  5  Mbar.  Cold  data  [1]  (filled  symbols)  and  Hugoniot  data  [19]  (unfilled 
symbols)  included. 
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WRIGHT  PATTERSON  AFB  OH  45433 

1  US  DEPT  OF  ENERGY 

IDAHO  OPERATIONS  OFHCE 
ATTN  RONALD  H  JOHNSON 
850  ENERGY  DR 
IDAHO  FALLS  ID  83401-1563 
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44  DIRECTOR  3 

SANDIA  NATL  LABS 
ATTN  E  S  HERTEL  JR  MS  0819 
LNKMETYK 
EHBARSISMS031 
ERIC  W  REECE  MS  0307 
DANIEL  P  KELLY  MS  0307 
LWEIRICKMS0327 

RTACHAUMS0425  33 

DLONGCOPEMS0439 

R REEDER 

D  HAYES  MS  0457 

JASAYMS0458 

WTEDESCHIMS0482 

J  SCHULZE  MS  0483 

J  SOUTHWARD 

PAUL  A  LONGMIRE  MS  0560 

TTRUCANOMS0819 

J  MICHAEL  MCGLAUN  MS  0819 

R  BRANNON  MS  0820 

JANG  MS  0821 

M  BOSLOUGH  MS  0821 

L  CHHABELDAS  MS  0821 

D  CRAWFORD  MS  0821 

JDICK 

M  FURNISH  MS  0821 
C  HALL  MS  0821 
W  REINHART  MS  0821 
P  STANTON  MS  0821 
P  TAYLOR  ORG  1432 
DKERNANORG1433 
C KONRAD 
KLANG 
M  KIPP  DIV  1533 
A  ROBINSON  DIV  1533 
P  YARRINGTON  DIV  1533 
JMCGLAWA  DIV  1541 
M  FORRESTAL  DIV  1551 
R  GRAHAM  DIV  1551 
R  LAFARGE  DIV  1551 
CHILLS  DIV  1822 
W  J  ANDRZEJEWSKI  DIV  2512 
D  MARCH!  DIV  2512 
W  VANDERMOLEN  ORG  2653 
B  LEVIN  ORG  7816 
R  O  NELLUMS  DIV  9122 
PO  BOX  5800 

ALBUQUERQUE  NM  87185-0307 


DIRECTOR 

LLNL 

ATTN  JOSEPH  V  REPA  MS  A133 
MICHAEL  O  SCHNICK  MS  F607 
EDWARD  J  CHAPYAK  MS  F664 
PO  BOX  1663 
LOS  ALAMOS  NM  87545 

DIRECTOR 

LLNL 

ATTN  RE  TIPTON  L35 

DBAUML35 

TMCABEEMS35 

MJMURPHY 

R  PIERCE  L122 

RROSINKYL122 

OJ  ALFORD  L122 

D  STEWART  L122 

TVIDLAKL122 

BR  BOWMAN  L122 

WDIXONL122 

ACMrrCHH.L 

JAMORIARTY 

RAHEINLE 

NC  HOLMES 

RPERRETL163 

WSHOTTSL163 

H  KRUGER  L178 

GPOMYKALL178 

M  SHANNON 

G  SIMONSON  L180 

A  SPERO  L180 

WTAOL282 

PURTIEWL282 

AHOLTL290 

JEREAUGHL290 

WJ  NELLIS  L299 

DWOODL352 

D  GAVEL  L495 

J  HUNTER  L495 

E  JOHANSSON  L495 

RMKUKLOL874 

GW  REPP 

PO  BOX  808 

LIVERMORE  CA  94550 
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8  DIRECTOR 
LLNL 

ATTN  DOUGLAS  R  FAUX  L125 
NORMAN  W  KLINO  L125H 
ROBERT  BARKER  L159 
MILTON  FINGER  L163 
MICHAEL  GERASSIMENKO  L178 
FRANK  A  HANDLER  LI  82 
STEPHEN  G  COCHRAN  L389 
BMDO  ROBERT  M  HALL 
PO  BOX  808 

LIVERMORE  CA  94551-0808 

34  DIRECTOR 
LANL 

ATTN  M  LUCERO  MS  A105 

DMANDELL 

KMARK 

J  P  RITCHIE  MS  B214  T14 
RDINGUSMSB218 
N  KRIKORIAN  MS  B228 
R  KIRKPATRICK  MS  B229 
R  THURSTON  MS  B229 
C  T  KLINGNER  MS  B294 
R  MILLER  MS  B294 
SJMOSSO 
B  SHAFER  MS  C931 
GGISLERMSD436 
C  RAGAN  MS  D449 
B  LAUBSCHER  MS  D460 
R  WELLS  MS  F607 
RKOPPMSF645 
R  STELLINGWERF  MS  F645 
C  WINGATE  MS  F645 
T  ADAMS  MS  F663 
R  GODWIN  MS  F663 
K  JACOBY  MS  F663 
W  SPARKS  MS  F663 
JSHANERMSF670 
G  CANAVAN  MS  F675 
R  GREINER  MS  G740 
J  HILLS  MS  G770 
B  HOGAN  MS  G770 
JBOLSTADMSG787 
J  WALSH  MS  G787 
R  DAVIDSON  MS  K557 
R  HENNINGER  MS  K557  N6 
TROLLETMSK574 
WDEALMSP915 
PO  BOX  1663 
LOS  ALAMOS  NM  87545 


12  DIRECTOR 
LANL 

ATTN  P  HOWE  MS  P915 
J  KENNEDY  MS  P915 
AROACHMSP915 
L  SCHWALBE 
W  HEMSING  MS  P940 
RKARPPMSP940 
E  POGUE  MS  P940 
J  MCAFEE  MS  P950 
D  PAISLEY  MS  P950 
L  PICKLESIMER  MS  P950 
RWARNESMSP950 
S  SHEFFIELD  MS  P952 
PO  BOX  1663 
LOS  ALAMOS  NM  87545 

3  ENERGETIC  MATERIALS  RSCHCTRTDOE 
NEW  MEXICO  INST  OF  MINING  &  TECH 
ATTN  DAVID  J  CHAVEZ 
LARRY  LIBERSKY 
FRED  SANDSTROM 
CAMPUS  STATION 
SOCORRO  NM  87801 

1  NASA 

LEWIS  RESEARCH  CENTER 
ATTNJFERRANTE 
CLEVELAND  OH  44135 

3  NASA 

JOHNSON  SPACE  CENTER 
ATTN  ERIC  CHRISTIANSEN 
JEANNE  LEE  CREWS 
FREDRICHHORZ 
MAIL  CODE  SN3 
2101  NASA  RDl 
HOUSTON  TX  77058 

1  APPLIED  RESEARCH  LAB 

ATTN  JEFFREY  A  COOK 
10000  BURNETT  ROAD 
AUSTIN  TX  78758 

1  GM  RESEARCH  LABS 
ATTNJRSMITH 
WARREN  MI  48090 
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5  ET  PROPULSION  LABORATORY 

ATTN  ZDENEK  SEKANINA 
PAULWEISSMAN 
BOB  WEST 
JAMES  ZWISSLER 

MARC  ADAMS  IMPACT  PHYSICS  GROUP 
4800  OAK  GROVE  DR 
PASADENA  CA  91 109 

1  MAXWELL  LABS 
S  CUBED  DIVISION 
ATTN  GERALD  A  GURTMAN 
PO  BOX  1620 
LA  JOLLA  CA  92037 

1  AMES  LABORATORY 

IOWA  STATE  UNIVERSITY 
ATTN  JH  ROSE 
AMES  lA  50011 

1  CALEORNIA  ST  UNIV 
NORTHRIDGE 
ATTN  JON  H  SHIVELY 
10343  4TH  ST  NW 
ALBUQUERQUE  NM  871 14 

2  CALTECH 

ATTN  ANDREW  P INGERSOIX  MS  170  25 
THOMAS  J  AHRENS  MS  252  21 
1201  E  CALEORNIA  BLVD 
PASADENA  CA  91125 

1  CALTECH 

ATTN  GLENN  ORTON  MS  169  237 
4800  OAK  GROVE  DR 
PASADENA  CA  91007 

1  DREXEL  UNIVERSITY 

ATTN  PHYSICS  DEPT 
32ND  &  CHESTNUT  ST 
PHILADELPHIA  PA  19104 

5  JOHNS  HOPKINS  UNTV 
APPLIED  PHYSICS  LAB 
ATTN  TERRY  R  BETZER 
ALVIN  R  EATON 
RICHARD  H  KEITH 
DALE  K  PACE 
ROGER  L  WEST 
JOHNS  HOPKINS  ROAD 
LAURH.MD  20723 


1  LOUISIANA  STATE  UNIVERSITY 

ATTN  ROBERT  W  COURIER 
948  WYLEDR 
BATON  ROUGE  LA  70808 

1  MIT  DEPT  OF  EARTH  ATMOS  AND 

PLANETARY  SCIENCES 
ATTN  HEIDI  B  HAMMELL  54  316 
CAMBRIDGE  MA  02139 

1  NC  STATE  UNIVERSITY 

ATTN  YASUYUKI HORIE 
RALEIGH  NC  27695-7908 

1  PENNSYLVANIA  STATE  UNIVERSITY 
ATTN  PHYSICS  DEPT 
UNIVERSITY  PARK  PA  16802 

5  SOUTHWEST  RESEARCH  INSTTIUTE 
ATTN  C  ANDERSON 
D  LIITLEFIELD 
SAMULLIN 
BCOUR  PALAIS 
J  WALKER 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78284 

1  STANFORD  UNIVERSITY 

ATTN  BRIAN  PETERSEN 
RAINS  16L 
704  CAMPUS  DR 
STANFORD  CA  94305 

1  TEXAS  A&M  UNIVERSITY 

PHYSICS  DEPARTMENT 
ATTN  DAN  BRUTON 
COLLEGE  STATION  TX  77843-4242 

1  UC  DAVIS 

ESrST  OF  THEORETICAL  DYNAMICS 
ATTN  EG  PUCKETT 
DAVIS  CA  95616 

1  UC  SANTA  BARBARA 

INSTITUTE  FOR  THEORETICAL  PHYSICS 

ATTN  F  GUINEA 

SANTA  BARBARA  CA  93 106 
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2  UNTV  OF  ALA  HUNTSVILLE 

AEROPHYSICS  RSCH  CTR 
ATTN  GARY  HOUGH 
DAVID  JLIQUORNIK 
PO  BOX  999 
HUNTSVILLE  AL  35899 

1  UNIV  OF  ALA  HUNTSVILLE 

CIVIL  ENGRNG  DEPT 
ATTN  WILLIAM  P  SCHONBERG 
HUNTSVILLE  AL  35899 

1  UNTV  OF  ALABAMA 

ATTN  A  EUGENE  CARDEN 
86  EASTERN  HILLS 
COTTONDALE  AL  35453 

1  UNIVERSITY  OF  ARIZONA 

ATTN  DAVID  H  LEVY 
TUCSON  AZ  85721 

1  UNIVERSITY  OF  CA  BERKLEY 
MECHNL  ENGRNING  DEPT 
GRADUATE  OFFICE 
ATTNKEZHUNLI 
BERKELEY  CA  94720 

1  UNIVERSITY  OF  CHICAGO 

DEPT  OF  THE  GEOPHYSICAL  SCIENCES 
ATTN  GH  MILLER 
5734SELUS  AVE 
CHICAGO  IL  60637 

1  UNIVERSITY  OF  COLORADO 
ATTNTMACLAY 
CAMPUS  BOX  431  NNT  341 
BOULDER  CO  80309 

4  UNIVERSITY  OF  DELAWARE 

DEPT  OF  MECHANICAL  ENGINEERING 

ATTN  PROF  J  GILLESPIE 

DEAN  RB  PIPES 

PROF  J  VINSON 

PROF  D  WILKINS 

NEWARK  DE  19716 

1  UNIVERSITY  OF  MARYLAND 

ATTN  PHYSICS  DEPT  (BLDG  082) 
COLLEGE  PARK  MD  20742 


1  UNIVERSITY  OF  PUERTO  RICO 
DEPT  CHEMICAL  ENGINEERING 
ATTN  LA  ESTEVEZ 
MAYAGUEZ  PR  00681-5000 

1  UNIVERSITY  OF  TEXAS 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN  ERIC  P  FAHRENTHOLD 
AUSTIN  TX  78712 

1  ADVANCED  TECHNOLOGY  INC 
ATTN  J  ADAMS 

PO  BOX  125 

DAHLGREN  VA  22448-0125 

2  AEROJET  ELECTRO  SYS  CO 
ATTN  WARHEAD  SYSTEMS 
JCARLEONE 

SKEY 
PO  BOX  296 
AZUSA  CA  91702 

2  AEROJET  ORDNANCE 
ATTN  P  WOLF 

G  PADGETT 

1100  BULLOCH  BLVD 

SOCORRO  NM  87801 

3  ALLIANT  TECHSYSTEMS  INC 
ATTN  T  HOLMQUIST  MNl  1  2720 
RSTRYK 

G  R  JOHNSON  MNl  1  2925 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  ALME  AND  ASSOCIATES 
ATTN  MARVIN  L  ALME 
6219  BRIGHT  PLUME 
COLUMBIA  MD  21044-3790 

1  APPLIED  RESEARCH  ASSOC  INC 
ATTN  JEROME  D  YATTEAU 
5941  S  MIDDLEFIELD  RD  SUITE  100 
LITTLETON  CO  80123 

1  APPLIED  RESEARCH  ASSOC  INC 
ATTN  FRANK  MAESTAS 
4300  SAN  MATEO  BLVD  SE 
ALBUQUERQUE  NM  87110 
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BATTELLE  1 

ATTN  ROBER  M  DUGAS 
7501  S  MEMORIAL  PKWY  SUITE  101 
HUNTSVILLE  AL  35802-2258 

BATTELLE 

ATTN  S  GOLASKI  1 

R  JAMESON 

TECHNICAL  LIBRARY 

505  JONG  AVE 

COLUMBUS  OH  43201 

2 

BOEING  AEROSPACE  CO 
SHOCK  PHYSICS  &  APPLIED  MATH 
ENGINEERING  TECHNOLOGY 
ATTNRHELZER 
T  MURRAY 
J SHRADER 

PO  BOX  3999  1 

SEATTLE  WA  98124 

BOEING  HOUSTON  SPACE  STN 

ATTN  RUSSELL  F  GRAVES  3 

BOX  58747 

HOUSTON  TX  77258 

BRIGS  CO 

ATTN  JOSEPH  E  BACKOFEN 
2668  PETERSBOROUGH  ST 

HERNDON  VA  2207 1-2443  3 

CALIFORNIA  RSCH  &  TECHNOLOGY 

ATTN  R  BROWN 

RFRANZEN 

DORPHAL 

51 17  JOHNSON  DR 

PLEASANTON  CA  94566  1 

CALIFORNIA  RSCH  &  TECHNOLOGY 

ATTNMMAJERUS 

PO  BOX  2229 

PRINCETON  NJ  08543  1 

COMPUTATIONAL  MECHANICS 

CONSULTANTS 

ATTNJAZUKAS 

PO  BOX  11314  1 

BALTIMORE  MD  21239-0314 


DEFENSE  TECHNOLOGY  INTERNATIONAL  INC 

ATTNDEAYER 

THE  STARK  HOUSE 

22  CONCORD  ST 

NASHUA  NH  03060 

DESKIN  RESEARCH  GROUP  INC 
ATTN  EDWARD  COLLINS 
2270AGNEWRD 
SANTA  CLARA  CA  95054 

El  DUPONT  DE  NEMOURS  &  CO 
ATTN  L  MINOR 
B SCOTT 

SECURITY  DIR  LEGAL  DEPARTMENT 
PO  BOX  1635 
WILMINGTON  DE  19899 

GE  DUVALL 

5814  NE82ND  COURT 

VANCOUVER  WA  98662-5944 

DYNAEASTCORP 
ATTN  PC  CHOU 
R  CICCARELLI 
WFLIS 

3620  HORIZON  DRIVE 
KING  OF  PRUSSIA  PA  19406 

DYNASEN 

ATTN  JACQUES  CHAREST 
MICHAEL  CHAREST 
MARTIN  LILLY 
20  ARNOLD  PL 
GOLETA  CA  93117 

ELORET  INSTITUTE 
ATTN  DAVID  W  BOGDANOFF  MS  230  2 
NASA  AMES  RESEARCH  CENTER 
MOFFETT  FIELD  CA  94035 

ENIG  ASSOCIATES  INC 
ATTNDJPASTINE 
13230 INGLESIDE  DR 
BELTSVILLE  MD  20705 

EXPLOSIVE  TECHNOLOGY 
ATTNMLKNAEBEL 
POBOXKK 
FAIRFIELD  CA  94533 
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1  GB  TECH  LOCKHEED 
ATTN  JAY  LAUGHMAN 
2200  SPACE  PARK  SUITE  400 
HOUSTON  TX  77258 

2  GB  TECH  LOCKHEED 

ATTN  LUCILLE  BORREGO  C23C 
JOE  FALCON  JR  C23C 
2400  NASA  ROAD  1 
HOUSTON  TX  77058 

2  GENERAL  RESEARCH  CORP 

ATTN  A  CHARTERS 
TMENNA 
PO  BOX  6770 

SANTA  BARBARA  CA  93160-6770 

2  GRC  INTERNATIONAL 

ATTN  TIMOTHY  M  CUNNINGHAM 
WILLIAM  M  ISBELL 
5383  HOLLISTER  AVE 
SANTA  BARBARA  CA  93 1 1 1 

5  INST  OF  ADVANCED  TECHNOLOGY 
UNIVERSITY  OF  TX  AUSTIN 
ATTN  STEPHEN  J  BLESS 
JAMES  CAZAMIAS 
HARRY  D  FAIR 
THOMAS  MKIEHNE 
NICK  LYNCH 
4030-2  WBRAKERLN 
AUSTIN  TX  78759 

1  KAMAN  SCIENCES  CORP 
ATTN  DENNIS  L  JONES 
2560  HUNTINGTON  AVE  SUITE  200 
ALEXANDRIA  VA  22303 

8  KAMAN  SCIENCES  CORP 
ATTN  J  ELDER 
RICHARD  P  HENDERSON 
DAVID  A  PYLES 
FRANK  R  SAVAGE 
JAMES  A  SUMMERS 
JAMES  S  WILBECK 
TIMOTHY  W  MOORE 
THYYEM 

600  BLVDS  SUITE  208 
HUNTSVILLE  AL  35802 


3  KAMAN  SCENCES  CORP 
ATTN  SHELDON  JONES 
GARY  L  PADEREWSKI 
ROBERT  GPONZINI 

1500  GRDN  OF  THE  GODS  RD 
COLORADO  SPRINGS  CO  80907 

4  KAMAN  SCIENCES  CORP 
ATTNNASITARI 
STEVE  R  DIEHL 
WILLIAM  DOANE 
VERNON  M  SMITH 

PO  BOX  7463 

COLORADO  SPRINGS  CO  80933-7463 

1  D  R  KENNEDY  &  ASSOC  INC 

ATTN  D  KENNEDY 
PO  BOX  4003 

MOUNTAIN  VIEW  CA  94040 

1  KERLEY  PUBLISHING  SERVICES 
ATTNGIKERLEY 

PO  BOX  13835 

ALBUQUERQUE  NM  87192-3835 

2  KTECH  CORPORATION 
ATTN  FRANK  W  DAVIES 
LARRY  M  LEE 

901  PENNSYLVANIA  NE 
ALBUQUERQUE  NM  87110 

1  LIVERMORE  SOFTWARE  TECH  CORP 

ATTNJOHALLQUIST 
2876  WAVERLY  WAY 
LIVERMORE  CA  94550 

1  LOCKHEED  ENGRNG  &  SPC  SCIENCES 

ATTN  E  CYKOWSKI  MS  B22 
2400  NASA  RDl 
HOUSTON  TX  77058 

1  LOCKHEED  MARTIN  MISSLE  &  SPACE 
ATTN  WILLIAM  R  EBERLE 
PO  BOX  070017 
HUNTSVILLE  AL  35807 
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3  LOCKHEED  MARTIN  MISSILE  &  SPACE 
ATTN  M  A  LEVIN  ORG  8 1  06  BLDG  598 
MRMCHENRY 

T  A  NGO  ORG  81  10  BLDG  157 
111  LOCKHEED  WAY 
SUNNYVALE  CA  94088 

4  LOCKHEED  MISSILE  &  SPACE  CO 
ATTN  JOHN  R  ANDERSON 
WILLIAM  C  KNUDSON 

S  KUSUMI 0  8111  BLDG  157 
J  PHILLIPS  0  54  50 
PO  BOX  3504 
SUNNYVALE  CA  94088 

1  LOCKHEED  MISSILE  &  SPACE  CO 

ATTN  R  HOFFMAN 
SANTA  CRUZ  FACILITY 
EMPIRE  GRADE  RD 
SANTA  CRUZ  CA  95060 

1  LOCKHEED  NASA  JSC 

SPACE  SCIENCE  BRANCH 
ATTN  JAMES  HYDE 
BOX  58561  MCB22 
HOUSTON  TX  77258 

1  LOCKHEED  MARTIN  AEROSPACE 
ATTNDRBRAGG 
PO  BOX  5837  MP  109 
ORLANDO  FL  32855 

1  MCDONNELL  DOUGLAS 
ASTRONAUTICS  CO 
ATTN  BL  COOPER 
5301  BOLSAAVE 
HUNTINGTON  BEACH  CA  92647 

1  ORLANDO  TECHNOLOGY  INC 
ATTN  DANIEL  A  MATUSKA 
PO  BOX  855 
SHALIMARFL  32579 

1  PHYSICAL  SCIENCES  INC 
ATTN  PETER  NEBOLSINE 
20  NEW  ENGLAND  BUS  CTR 
ANDOWERMA  01810 


3  PHYSICS  INTERNATIONAL 
ATTNRFUNSTON 
G  FRAZIER 
L  GARNETT 
PO  BOX  5010 
SAN  LEANDRO  CA  94577 

1  ROCKWELL  INTERNATIONAL 

ROCKETDYNE  DIVISION 
ATTN  J  MOLDENHAUER 
6633  CANOGA  AVE  HB  23 
CANOGA  PARK  CA  91303 

1  ROCKWELL  INTERNATIONAL 
ROCKETDYNE  DIVISION 
ATTNHLEIFER 
16557  PARK  LN  CIRCLE 
LOS  ANGELES  CA  90049 

1  ROCKWELL  MISSILE  SYS  DIV 

ATTNTNEUHART 
1800  SATELLITE  BLVD 
DULUTH  GA  30136 

1  S  CUBED 

ATTN  R  SEDGWICK 
PO  BOX  1620 
LA  JOLLA  CA  92038-1620 

1  SAIC 

ATTN  GREGORY  J  STRAUCH  MS  264 

1710GOODRIDGEDR 

MCLEAN  VA  22102 

1  SAIC 

ATTN  MICHAEL  W  MCKAY 
10260  CAMPUS  POINT  DR 
SAN  DIEGO  CA  92121 

1  SHOCK  TRANSIENTS  INC 
ATTN  DAVID  DAVISON 
BOX  5357 

HOPKINS  MN  55343 

2  SIMULATION  &  ENG  CO  INC 
ATTN  ELSA  I  MULLINS 
STEVEN  EMUULINS 

8840  HWY  20  SUITE  200  N 
MADISON  AL  35758 
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SOUTHERN  RESEARCH  INSTITUTE 
ATTN  LINDSEY  A  DECKARD 
DONALD  PSEGERS 
PO  BOX  55305 

BIRMINGHAM  AL  35255-5305 

SRI  INTERNATIONAL 
ATTN  JAMES  D  COLTON 
MOHSENSANAI 
DR  L  SEAMAN 
DONALD  A  SHOCKLEY 
333  RAVENSWOOD  AVE 
MENLO  PARK  CA  94025 

TELEDYNE  BROWN  ENGR 
ATTN  JIM  W  BOOTH 
MARTIN  B  RICHARDSON 
SHERRY  RIDDLE 
PO  BOX  070007  MS  50 
HUNTSVUXE  AL  35807-7007 

ZERNOW  TECHNICAL  SVCS  INC 
ATTN  LOUIS  ZERNOW 
425  W  BONITA  AVE  SUITE  208 
SAN  DIMAS  CA  91773 


ABERDEEN  PROVING  GROUND 

37  DIR,  USARL 

ATTN:  AMSRL-WM,  I.  MAY 

AMSRL-WM-T,  W.  F.  MORRISON 
AMSRL-WM-TA, 

W.  BRUCHEY 
J.DEHN 

G.  FILBEY 
W.  GILLICH 
W.  A.  GOOCH 

H.  W.  MEYER 

E.  J.RAPACKI 
AMSRL-WM-TB, 

R.FREY 

J.  STARKENBERG 
AMSRL-WM-TC, 

W.S.DEROSSET 
T.  W.  BJERKE 

R.  COATES 

F.  GRACE 

K.  KIMSEY 
M.  LAMPSON 
D.  SCHEFFLER 

S.  SCHRAML 
B.  SORENSEN 

R.  SUMMERS 
W.  WALTERS 

AMSRL-WM-TD, 

A.  M.  DIETRICH 
K.  FRANK 
J.  HARRISON 
M.  RAFTENBERG 

G.  RANDERS-PEHRSON 
J.  SANTIAGO 

M.  SCHEIDLER 

S.  SEGLETES  (3  CP) 

J.  WALTER 

T.  WRIGHT 

AMSRL-WM-PB,  A.  ZIELINSKI 
AMSRL-WM-PD,  G.  GAZONAS 
AMSRL-WM-WD,  A.  PRAKASH 
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1  ABTEILUNGFUERPHYSIKALISCHE 
CHEMIE 

MONT  ANUNTVERSITAET 
ATTN  E  KOENIGSBERGER 
A8700LOEBEN 
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